Casimir stress on lossy magnetodielectric spheres 



Christian Raabe, Ludwig KnoU, and Dirk-Gunnar Welsch 
Theoretisch-Physikalisches Institut, Friedrich-Schiller-Umversitdt Jena, Max-Wien-Platz 1, D-07743 Jena, Germany 

(Dated: February 1, 2008) 

An expression for the Casimir stress on arbitrary dispersive and lossy linear magnetodielectric 
matter at finite temperature, including left-handed material, is derived and applied to spherical 
systems. To cast the relevant part of the scattering Green tensor for a general magnetodielectric 
sphere in a convenient form, classical Mie scattering is reformulated. 

PACS numbers: 



I. INTRODUCTION 

Spheres and spherical shells display enough symme- 
try to allow for a quite explicit solution of the Casimir 
problem (for the earliest calculation, see 0). Especially 
two intriguing but nowadays obsolete ideas have substan- 
tially stimulated the study of spherically symmetric se- 
tups. Firstly, the idea that the Casimir effect might be 
responsible for the postulated Poincare stress that mirac- 
ulously holds the charge distribution of a classical elec- 
tron together ^2], and secondly the idea that it might be 
related to sonoluminescence 0,1^, i.e., the emission 
of light flashes by small air bubbles in water under the 
action of ultrasonic waves. 

Unfortunately, in the available work on the topic 
(see, e.g., 0, H H ^> 01) lossless materials are as- 
sumed, dispersion is ignored, or the influence of matter 
is summarized in perfect conductor boundary conditions 
P, 0, 0. Needless to say that such approximations 
should be expected to distort the physical picture more 
or less drastically, not to mention their tendency to cause 
'unreal' divergence problems. In what follows, we al- 
low for dispersing and absorbing (causal) sphere material 
characterized by both a complex frequency-dependent 
permittivity and a complex frequency-dependent perme- 
ability, thus including in the calculations also left-handed 
materials, which have been of increasing interest. To 
our knowledge, Casimir forces on magnetodielectric bod- 
ies have been studied only for non-dispersing and non- 
absorbing material 0, commonly subject to the condi- 
tion that the product of the permittivity and the perme- 
ability is also uniform in space B, IS Uo': 111 . By using the 
quantization scheme given in Ref. |l4] in our calculations, 
we first extend the recently derived basic formula 15] for 
the Casimir force for causal dielectric matter to causal 
magnetodielectric matter. We then apply the theory to 
spherical structures. 



II. QUANTIZATION SCHEME 

Within the mentioned quantization scheme, the macro- 
scopic (medium-assisted) electromagnetic field operators 
are all expressed in terms of suitable bosonic basic fields 
via the classical Green tensor, which satisfies the equa- 



tion 

^2 

Vx K(r, V X G(r, r', oj) -£(r, uj) G(r, r', uj)^S(r, r') 

, (1) 

(e ~ complex permittivity, k = fjL ^ - complex recipro- 
cal permeability) and the boundary condition at infinity. 
In the absence of additional charges and currents, the 
positive-frequency components of the electric field oper- 
ator are then given by 

E(r, uj) = ^^io^ j cPt' G(r, r', c^)j^(r, uj), (2) 

from which the other field operators can be derived by 
using Faraday's law 

B(r, Lo) = {iuj)-^V X E(r, uj) (3) 

and the constitutive relations [kq — ^Iq^] 

D(r,c^) = £o£(r,w)E(r,c^)-f PN(r,c^), (4) 
H(r,a;) = kok{v,uj)'B{v,uj) ~M^{r,uj). 

In the above, the frequency-domain Langevin noise quan- 
tities (carrying an index N) are connected with the funda- 
mental bosonic fields f\{Y,Ljj) {X — e,m) by the relations 

jj^(r,w) = -iujP^{r,uj) + V X M^{r,^), (5) 
P^{r, uj) = i [feoIme(r,w)/7r]^^^fe(r,w), (6) 

^ 1 /2 ^ 

MN(r,w) = [-hKolinK{r,uj)/TT] fm(r,w). (7) 

As a consequence, the electromagnetic field operators sat- 
isfy the usual commutation relations of QED. For exam- 
ple, the electric field operator (in the Schrodinger picture) 
reads as 

poo 

E(r) = / dwE(r,tj) +H.C., (8) 
Jo 

with E(r, w) from Eq. ^ together with Eqs. 10-0. 
Writing down the corresponding formulas for the other 
field operators is straightforward. The consistency of the 
method relies heavily on the facts that the Green tensor, 
which is a response function like £(r, w) and «;(r, w), is 
holomorphic in the upper lu half-plane, obeys the reci- 
procity relation 

G{r,r',u;)^ G^{r',r,u;) (9) 



2 



(the superscript T denotes matrix transposition) , and the 
integral relation 



d^s I [G(r, s, w) X V,] Ini k(s, w) [V, x G* (s, r', w)] 

H — Y G(r, s, Li;)Im£(s, lu) G*(s, r', a;) > =Ini G(r, r', w) 

(10) 

(for details, see ^3)- Note that the convention 

G(r, r', w) X V' = - V' x G'^(r, r', oj) (11) 
has been used. 

III. CASIMIR STRESS TENSOR 

For a given quantum state, the stress tensor, i.e., the 
(yet unrenormalized) Casimir stress, can be obtained 
from the correlation function 



T(r, r') = (D(r) ® E(r')) + (B(r) ® H(r')) 

(D(r)®E(r')) + (B(r)®H(r')) 



ilTr 



(12) 



in the coincidence limit r' r. In carrying out this limit, 
we have to drop - in accordance with the Casimir effect's 
very definition - the bulk part of the Green tensor, which 
is singular but independent of geometry p^ . to get a 
physical (finite) value of the Casimir force per unit area. 
Note that the resulting force formula rigorously applies 
only to points outside matter, which is the case under 
consideration. 

Let us calculate the thermal-equilibrium Casimir force 
at temperature T. Recalling the bosonic character of the 
fundamental fields {x{r,uj) and assuming them to be ex- 
cited in thermal states, one quickly finds, in close analogy 
to Ref. fii, that 



Sxx'S{u;~Lu')Sir,r'), (13) 



coth — - 1 

\2kBTj 



Sxx'Siu;-Lu')Sir,r'), (14) 



{ix{r,u;)^h'{r',uj')) = {il{r,u;)<g>ix'{r' ,u;')) = 0. (15) 

By using Eq. ^ for the electric field and the related 
equations for the other fields [together with Eqs. ©-Q] 
and Eqs. (|13|I - H15(I . and by employing the Green tensor 
properties Q and (|10|l . it follows by a cumbersome but 



straightforward calculation that Eq. (|12(l may be rewrit- 
ten as 



r(r,r') = 0(r,r')-iJTr6>(r,r'), (16) 



where 



X Im< ^e(r, lu) G(r, r', lu) 



Vx G(r,r',a;) x V nir', lu) }. (17) 



Note that this result differs - apart from the differ- 
ent analytical form of the Green tensor - from the 
'non-magnetic' result in the permeability that ap- 
pears in the second term in the curly brackets. By re- 
placing the full Green tensor with its scattering part, 
G(r,r',w)i— > G'^'^^'(r, r', w), and taking the coincidence 
limit r' r, Eq. l(TH)l together with Eq. l(T7|l is the sought- 
after basic formula for the Casimir force per unit area 
for an arbitrary arrangement of lossy magnetodielectric 
bodies. With the bulk Green tensor removed and appli- 
cation of the standard rule lu i— > lu+iO to handle free-space 
regions, no divergence problems arise in performing the 
integral in Eq. ((T7|l . 



IV. SPHERE GREEN TENSOR 

To apply the theory to spheres, we need a convenient 
expression for the sphere scattering Green tensor - a 
purely classical problem, which was first treated by Mie 

eand is now standard textbook material (see, e.g., 
01). Here, we give a formulation that is valid for 
arbitrary complex permittivities and permeabilities that 
vary radially in a stepwise fashion. (Note that ampli- 
fying materials are in fact also allowed.) The required 
techniques can be gathered from the mentioned volumes, 
but see also 20^, JlJ . According to the method of De- 
bye's potentials, transverse solutions of the vector wave 
equation 



V X V X F(r) - fc2F(r) = 



(18) 



are given, for each Debye potential u that solves the 
scalar Helmholtz equation (A -I- k'^)u = 0, by vector func- 
tions F = V X ru and F = V x V x ru. In addition, all 
longitudinal vector functions F — Vu solve Eq. (|18|l for 
k — 0. Note that the free-space eigenfunctions of the 
(semi-bounded) operator Vx Vx obtained from Eq. l(TH|) 
under the restriction that they be bounded everywhere 
(including infinity) have real k and are complete |2l| . In 
free space, an orthonormalized and complete set of scalar 
(continuum) eigenfunctions u reads 



ui^ir;k)^\l-kji{kr)Yimie,(b) (fc > 0), (19) 
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with 1 = 0,1,2 .. ., m =—l, . . . ,1 (ji - spherical Bessel func- 
tion, - spherical harmonic |23|). Employing them as 
Debye potentials and introducing appropriate normaliza- 
tion factors yields the orthonormalized vector eigenfunc- 
tions as 

L;™(r; k) = ■ivu;„(r; k) 



= \l-[kfiikr)rYU0,cP)/r+ji{kr)VYi^{e,^)],{2O) 



where the primed sum begins with 1 = 1. Note that the 
M and N functions represent, respectively, TE and TM 
(to r) partial waves, and ru plays the role of a Hertz 
vector. 

To find the sphere scattering Green tensor, we will 
make use of the vector functions just constructed. For 
this purpose, let us first consider wave propagation in 
a homogeneous, isotropic bulk material, so that Eq. 
simplifies to 



Mimir;k) = [l{l + l)] 5Vxru/™(r;fc) 



ttI{1 + 1) 



kji{kr)r x VY, 



(21) 



V X V X G('^""^)(r,r',w) 

- k\u) G('^"i'^)(r, r', Lo) = fi{cu)S{r, r'), (25) 

with arbitrary complex 



Ni™(r;fc) = - V X Mi™(r;fc) 

k 



^ ^1^^^^^ {l{l + mkr)rYUO,<t>)/r^ 

+ [krjiikr) + jiikr)] VYi„,{e,c^)} , (22) 

where for the M and N functions the value / = has to be 
excluded from consideration. Throughout the literature 
it is invariably assumed that the Debye potential method 
does not miss any eigenf unctions. We also do so and 
assume the completeness relations 

/•OO 

^/ dfcW(r;fc)®LL(r';fc) = Jll(r,r'), (23) 
Im ^0 

J2 / & [Mi™(r;fc)0ML(r';fc) 

+N,„,(r; k) ® N*„(r'; k)] = S^{r, r'), (24) 



(26) 



The bulk-material Green tensor can then be expressed 
straightforwardly in terms of the vector functions H21|) 
and ^ according to 



G'(^"^^)(r,r',c.)M-(c.) = -^^ 

+ J2 I TT^^r^rr [M,™(r;fc')®ML(r';fc') 



Im 



k'^-k^uj) 



+N,™(r;fc')®Nr„(r';fc')], (27) 



and the integral can be evaluated by means of the residue 
theorem (with due care of 'static' poles in the N (g) N* 
terms 17]) to obtain [k^kiui)] 



(r — r) 



^2 ^(r, r') + ^ E' { \^irn{v; k) ® ML(r'; k) + Ni™(r; k) ® N.'Jr'; k) 

Im 

-f [M/™(r; k) ® M'„(r'; k) + N;„(r; fc) ® N'„(r'; fc)] e{r - r')} . (28) 



Here, the 'bullet' symbol (*) denotes complex conjuga- 
tion of the fc-independent factors only, and the tilde 
symbol ("') means that in the definitions 1)21(1 and 122(1 
the spherical Bessel functions ji have to be replaced 
with the 'outgoing' and 'incoming' spherical Hankel func- 
tions hi ~ hf'^ and hi = hf'^ for Im /c > and Im /c < 0, 
respectively, to ensure amplitude decay as |r — r'j^cx). 
The advantage of this notation is that the sign in 



k = -iiLo^feJl/c can be chosen freely. In this way, the for- 
mulas obtained for the Green tensor apply also to left- 
handed and amplifying material. 

The terms within the curly brackets in Eq. ((28(1 cor- 
respond to TE and TM partial waves that are solutions 
to the homogeneous version of the differential equation 
((25(1 with the dispersion relation ((26(1 . In a spherically 
layered medium (in contrast to the bulk material con- 
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sidered so far), additional waves that arise from reflec- 
tion at and transmission through layer interfaces must 
be taken into account. They just form the scattering 
part of the Green tensor that we are interested in. The 
total field must satisfy the well-known continuity con- 
ditions at the layer interfaces, the required amplitude 
matching can be done for each partial wave separately. In 
fact, only the solutions ~ kSi{kr) of the radial Helmholtz 
equation [^+l^-iil^+k^]Si{kr)= are involved in 
this matching process. To be more specific, Si{kr) must 
be a (bounded) superposition of the different kinds of 
spherical Bessel functions (for complex k) in every layer 
such that 

• for TE waves continuity of 

krSiikr) and fi'^k-^ [rSi{kr)] (29) 



and for TM waves continuity of 

d 

ekrSi(kr) and k—[rSi(kr)] (30) 
or 



is ensured at spherical interfaces r = const. 

The simplest case is a homogeneous sphere (index 1) 
of radius i? in a homogeneous environment (index 2). 
If both r and r' are outside the sphere and R< r <r' 
is valid, it can be concluded from the above mentioned 
conditions of continuity that the reflection coefficients 
related to the terms with 6{r' — r) in Eq. I|28|) are given 

by 



'J, 21 



and 



„TM 
W.21 



jijkiR) d[Rji(k2R]] _ j,{k2B.) a[R.j,{kiR)] 

^2 dR. fi-i OR 

jijkiR) d[Rhi{k2R)\ hi{k2R) dlRjiikiR)] 

fi2 dR Ml 9R 



ji(kiR) dlRji{k2R.)] _ ji{k2R) d[Rji{kiR)] 

£2 dR €i dR 

3i(kiR) d[Rki(k2R)] hi(k2R) d[Rji(kiR)] 

£2 dR £1 dR. 



(31) 



(32) 



for TE and TM waves, respectively. Introducing the def- 
inition 



V[jiikiR),hi{k2R);a 
= det 



ji{kiR) hi{k2R) 
a^^dji{kiK)/dR a^^dhi{k2R)/dR 



(33) 



we may write Eqs. (|31|l and (|32|l in the compact form of 
V[RJliklR),RJl{k2R)■,^i] 



and 



TE 
'^1,21 



TM 
' i,21 



V[Rji{kiR),Rhi{k2R);^i] 

V[RjiikiR),Rji{k2R);e] 
'v[Rji{kiR),Rhi{k2R);e] ' 



(34) 



(35) 



respectively. Note that the so-called Mie resonances of 
the sphere can be found by studying these expressions in 
the complex uj plane. The contributions to the scattering 
Green tensor for the case r>r'>R [related to the terms 
with 9{r — r') in Eq. (|28|l ] need not be figured out sep- 
arately, but may also be found from the condition that 
the scattering Green tensor has to be continuously differ- 
entiable at r = r' (since the bulk part alone accounts for 
the necessary singularity). The scattering Green tensor 
can therefore be given, for the case that both arguments 
r and r' are outside the sphere, in the form of 

(r, r^ ^) M,„.(r; fc2) ® ML(r'; fc2) 



Im 



+ r™Ni™(r;fc2)®Nr„(r';fc2)}. (36) 

By construction, G^^^-* (r, r', w) is continuously differen- 
tiable at r = r', vanishes for |r| — > oo and/or |r'| — ^ oo, and 
satisfies the condition of reciprocity. Note that in fact 
' = , because of the addition theorem 

' 2/ + 1 

Yi^i0,<p)YC^{e',.p') ^ ^ Piicosi;) gR (37) 

m—~l 

of the spherical harmonics {ip is the angle between the 
primed and unprimed directions) . It is worth noting that 
the form of Eq. H36|l does not change when the sphere is 
not homogeneous but consists of radial layers. Clearly, 
the reflection coefficients (|34|l and (|35|) must then be re- 
placed by generalized ones that take 'subsurface' reflec- 
tions into account. These can be calculated recursively 
by standard methods the basis of the matching 

conditions and 



V. CASIMIR STRESS ON A SPHERE 

With Eq. H3fci|l at hand, we can now evaluate the 
stress tensor according to Eqs. (I16|l and H17|) [we set 
G gr*5cat ^ q'-^^'i ^ to approach the sphere surface from 
outside], on using the symmetry between the M and N 
functions and assuming free space around the sphere, i.e., 
£2 = ^2 = 1, k2=ij/c + i0: 



T(r,r)= 6»(r,r)-iiTr0(r,r), 



(38) 



0(r,r) = A /li--coth(^) Re5:'(r™ +r™) 

M,„,(r; fc2)<»ML(i"; fc2) + Ni™(r; fc2)®Nf„(r; ^2) . 

(39) 

By symmetry, the relevant stress tensor element for the 
spherical setup is 



Ter = 6»e.r - Tr /2 



(40) 
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(e^ = v/r). By taking into account that the addition 
theorem (|37f) imphes the relation 



-—I 



21 + 1 

477 ' 



(41) 



and making use of Eqs. (|21|) and H22|) it is not difficult to 
see that the relations 



GrMi^ir- k) (g) ML(r; fc)e,. = 0, 



(42) 



(2/ 



fef(fcr) ?(/ + l) 
2^2 .2 ' 



(43) 



M,„(r;ft)ML(r;fc) = (2Z + l) 



^[(fcr) 
27r2 



fc^ (44) 



and 



ENi„.(r;fc)Nr^(r;fc) 

m 

2/ + 1 



l{l + \)ht{kr)+ krh'iikr) + hi{kr) W45) 



27r2r 



are valid. Substituting Eqs. lO-lllSl) in Eq. (123), we 
obtain from Eq. 14U|I the Casimir force per unit area on 
the surface of a sphere of radius R: 



RR 



auj ijj coth 

yiksT 



ReE(2^ + l)(r;^|i+r™)F,, (46) 



where 



Fi = 



^^^^^^ hf{k2R)~kfhUk2R) 



i?2 



i?2 



k2Rh'i{k2R) + hi{k2R) 



(47) 



Note that the frequency integral can be further processed 
by the standard contour deformation to the imaginary 
frequency axis, thereby leading to a sum over residues 
(Matsubara frequencies) when T > 0. 



VI. SUMMARY 



We have derived a general expression [Eqs. H16|l and 
(|17|l ] for the Casimir stress that is valid for arbitrary dis- 
persing and absorbing magnetodielectric material, and 
have evaluated it for (layered) spheres. In particular, 
Eq. H46|l [together with Eq. H47|l ] is applicable to spheres 
made of left-handed material, where interesting results 
can be expected when the relevant Mie resonances corre- 
spond to spectral regions in which the material behaves 
left-handed. Since the outlined Green tensor construc- 
tion also allows for linear amplification, the theory of 
Casimir forces may also be extended to (linearly) ampli- 
fying bodies. For this purpose, the calculations can again 
be based on the quantization scheme given in Ref . |l4| , if 
appropriate modifications (of the definitions of the fun- 
damental bosonic fields) are taken into account. Such 
modifications will effectively add to Eq. (|17|1 another in- 
tegral, over those frequency intervals and spatial regions 
which can contribute to the amplification. 
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